Exact ground state of the Shastry-Sutherland lattice with classical Heisenberg spins 
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An exact analytical solution of the ground state problem of the isotropic classical Heisenberg 
model on the Shastry-Sutherland lattice in external magnetic field H is found for arbitrary ratio of 
diagonal and edge exchange constants J2/ Ji . The phase diagram of this model in the ( J2/ Ji, H/ J\) 
plane is presented. It includes spin-flop, spin-flip and umbrella phases. The magnetization curves 
are found to be linear until saturation. It is shown numerically that the inclusion of the easy-axis 
anisotropy into the model leads to the appearance of the 1/3 magnetization plateau, corresponding 
to the collinear up-up-down spin structure. This explains the appearance of the 1/3 magnetization 
plateau in rare earth tetraborides RB4. In particular, magnetization curve of the compound H0B4 
is explained. 

PACS numbers: 75.10.Hk, 75.30.Kz, 75.40.Cx, 75.60.Ej 
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I. INTRODUCTION 

Shastry-Sutherland lattice (SSL) was first introduced 
in the work of Shastry and Sutherland 1 as a purely theo- 
retical example of a two-dimensional frustrated spin sys- 
tem. SSL is a square lattice with classical Heisenberg, 
quantum Heisenberg, or Ising spin at every lattice site 
i, with antiferromagnetic (AFM) exchange J\ along the 
edges, and AFM exchange J2 along certain diagonals, as 
shown in Fig. 1. Its Hamiltonian (in the presence of the 
external field H directed along the z-axis) is 



Ji 



edges 



diagonal 



and the Hamiltonian normalized by J\ is 

h^^=J2 Si s i + p E s ^ - h J2 ( 2 ) 



edges 



diagonal 



where we have introduced the notations 
J 2 , H 



(3) 



Surprisingly, ten years after the work of Shastry and 
Sutherland, SSL has been experimentally realized in 
SrCu2(BC>3)2 2 , which has a layered structure, with each 
Cu 2+ ion carrying spin S — 1/2. If exchange interactions 
with only two spheres of nearest neighbors are included, 
the magnetic lattice of SrCu2(B03)2 is equivalent to SSL. 
At low temperatures SrCu.2 (60*3)2 exhibits a sequence of 
magnetization plateaus at fractional values of the satura- 
tion magnetization M s 3-5 . A number of theories of this 
phenomenon has been proposed 6-8 . 

Similar fractional magnetization plateaus have been 
observed recently in rare earth tetraborides 9-12 RB4, 
where R=Tm, Er or Ho, where the rare earth ions also 
form layered structure equivalent to SSL. The important 
difference is that while Cu 2+ ions in SrCu2(B03)2 have 
spins s = 1/2, the rare earth ions in RB4 systems have 
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FIG. 1. (Color online) Shastry-Sutherland lattice. The red 
(gray) triangle marks the elementary cluster of the SSL. 



large spins, which can be treated as classical ones. These 
compounds also possess a strong easy-axis magnetocrys- 
talline anisotropy. 

The discovery of magnetization plateaus in RB4 com- 
pounds led to a number of theoretical and computa- 
tional studies of SSL with classical Heisenberg 13-15 and 
Ising 16-21 spins. Deep understanding of these two models 
is vital for explaining the peculiar magnetization curves 
of RB4 systems, as they can serve as the foundation 
stones on which more complicated models with additional 
interactions can be built. A major breakthrough for 
Ising SSL in external magnetic field came very recently 
as its exact ground state has been found analytically 19 . 
This model gives a single M/M s = 1/3 magnetization 
plateau, which corresponds to the so-called up-up-down 
(UUD) phase (Fig. 2, lower left). On the other hand, 
a Monte Carlo simulation for the classical Heisenberg 
SSL 13 found magnetization curves with no steps. How- 
ever, when the easy-axis anisotropy was included in the 



2 



model 14,15 , the 1/3 UUD plateau appeared for a cer- 
tain range of Jij J\- Magnetization steps other than 1/3 
do not appear for either Ising or Heisenberg SSL. Addi- 
tional exchange or dipolar interactions 16,19,21 or lattice 
disorder 14 were employed to account for those plateaus. 
If these additional interactions are included, Ising model 
on SSL essentially succeeds in explaining the appearance 
of fractional magnetization plateaus in RB4. A typical 
zero-temperature magnetization curve M(H) of the Ising 
SSL has the "staircase" shape, namely it consists of hor- 
izontal magnetization steps (including the M = and 
M = M s ones) separated by first-order phase transitions 
(vertical segments of the M(H) curve). This is very sim- 
ilar to the experimental M(H) curve of TmB4 (Ref. 11). 
Anisotropic Heisenberg SSL 14,15 , on the other hand, gives 
smoother M(H) curves with inclined regions, which were 
experimentally observed for H0B4 (Ref. 12) and possibly 
ErB 4 (Ref. 9). 

While the exact ground state of the Ising SSL in mag- 
netic field has been found, the complete understanding 
of the classical Heisenberg SSL is still lacking. In par- 
ticular, its exact ground state has not been determined, 
except for special cases H = (Ref. 1) and J2/J1 = 2 
(Ref. 13). The Monte Carlo simulations 13-15 mainly fo- 
cused on the special point J%/ J\ — 2 or its vicinity, while 
the phase diagram in the (J2/J1, H/Ji) plane has never 
been published. The present paper is an attempt to clar- 
ify these issues. Its goal is to study in detail the ground- 
state problem of the isotropic classical Heisenberg SSL 
and to determine the phase diagram of this model in the 
(p, h) plane. The very interesting problem of the clas- 
sical Heisenberg SSL with easy-axis anisotropy is also 
addressed briefly in the present paper. 

The paper is organized as follows. Section II intro- 
duces different possible magnetic structures of the SSL 
and presents its phase diagram in the (p, h) plane. In sec- 
tion III the exact ground state of the SSL is calculated 
and it is proven that the three phases of the previous 
section are indeed the ground state spin structures. Sec- 
tion IV checks the exact result with numerical simulation 
and examines the effect of the uniaxial anisotropy on the 
magnetization curves. It is followed by a conclusion. 

II. MAGNETIC STRUCTURES AND THE 
PHASE DIAGRAM 

In this and the following sections we consider the prob- 
lem of the ground state of the SSL as the function of two 
parameters: h = H/Ji and p = J-zj J\. In the absence 
of external magnetic field (h = 0) the problem has been 
solved in Ref. 1. For p < 1 (J2 < Ji) the ground state is 
the regular Neel AFM state (Fig. 2, upper left) with the 
energy per lattice site 

e Ned = -2+£. (4) 

The Neel state satisfies all J\ exchange interactions, but 
not J2 ones. For p > 1 (J2 > Ji) the ground state is 
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FIG. 2. (Color online) Neel, spin-flip, UUD and spiral (for 
J2/J1 = 2) spin structures. Black and white circles denote 
spins pointing up and down respectively. Red, green and blue 
circles denote three different directions of the spiral structure, 
with angle 120° between each two directions. The magnetic 
unit cells are shown with purple (gray) rectangles. 



the so-called spiral state. The angle between neighbor- 
ing spins is 7r — Acf) along the edges, and 2Acf> along the 
diagonals, where A<j> — cos~ 1 (l/p) = cos _1 ( J\j J2). The 
energy of this structure is 

1 P IK\ 
^Spiral — ~Z ■ [0) 

P ^ 

Such configuration can be constructed in an infinite num- 
ber of ways, leading to a degeneracy 1,13 . In general, it is 
incommensurate with the crystal lattice, but for chosen 
values of p, namely for cos _1 (l/p) = irm/n with integer 
to, n, periodic spirals can be realized. One possible spiral 
configuration for p = 2 is presented in Fig. 2, lower right. 

The special case p = 2 (J2/J1 = 2) for h > has been 
solved in Ref. 13. In this case the Hamiltonian (1) pos- 
sesses an additional degeneracy, and there is an infinite 
number of spin structures which share both total energy 
and the magnetization with the umbrella structure in- 
troduced below. The generic case h > is a bit more 
complicated. First we introduce several candidate spin 
structures which correspond to the local extrema of the 
total energy and present the phase diagram of the SSL in 
the (p, h) plane. In the next section we prove rigorously 
that the three structures considered (spin-flop, spin-flip 
and umbrella) are indeed the lowest-energy structures in 
the respective regions of the (p, h) plane. 

UUD structure, shown in Fig. 2, has energy and mag- 




FIG. 3. (Color online) Umbrella structure for p = 1.5, h — 
3 calculated numerically on the 6x6 lattice with periodic 
boundary conditions. The arrows show the xy components of 
spins. The red (gray) circles show the z components of spins 
(all equal for the umbrella structure) . 



netization: 

zuud = -\ - ~ - -, M UUD = l/3. (6) 
3 6 3 

Here and in the following we define M as the magneti- 
zation per lattice site, so that M s = 1. The spin-flop 
structure is the Neel structure with all spins tilted by 
the angle 9' = tt/2 — 9 = sin -1 (h/8) out of the xy plane. 
This solution exists for h < 8, and its energy and mag- 
netization are given by 

epiop = -2 + | - ^, Mpiop = cos(0) = ^. (7) 

The spin-flip (ferromagnetic) structure (Fig. 2, upper 
right) has all spins aligned along the magnetic field. It's 
energy is 

e F u P = 2 + | - h, M FHp = 1. (8) 

At h = 8 the spin-flop structure turns into the spin-flip 
structure in a continuous fashion. Since M(h) is contin- 
uous at the point h — 8, but dM(h)/dh is not, this is a 
second-order phase transition. 

Another possible structure is the umbrella structure, 
proposed in Ref. 13. It is essentially a tilted spiral struc- 
ture. In the umbrella structure the spherical angles 9i of 
all spins are equal, and the angles cj>i are distributed like 
in the spiral structure above. The energy for given 0, 
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FIG. 4. (Color online) Umbrella structure for p — 1.5, h — 3 
calculated numerically on the 12 x 12 latice with periodic 
boundary conditions. The arrows show the xy components of 
spins. The red (gray) circles show the z components of spins 
(all equal in this case). 



is 

e(6, Acf>) = 2 (- sin 2 9 cos(A< ? !)) + cos 2 9) + 

| (sin 2 9 cos(2 A0) + cos 2 9) - h cos 9. (9) 

Minimization with respect to 9, Acf> gives A<fi — 
cos^ 1 (l/p) (it does not depend on h), and 

M Umb = cos 9 = kp (10) 
2{p + l) 



f-Umb — 7T ; -J- (J-1J 

The umbrella structure exists for p > 1 and h < h max — 
2 (p + l) 2 / p, has energy lower than the spin-flop one, and 
turns into the spin-flop structure at p = 1. At h = h max 
it becomes the spin-flip structure. Both phase transitions 
are of the second order. In Fig. 3 and Fig. 4 numerically 
calculated umbrella structures for (p, h) — (1.5,3) and 
(1.2, 3) respectively are presented (See section IV below 
for details). 

The phase diagram of the classical isotropic SSL is 
shown in Fig. 5. The solid lines mark the lines of 
the second order phase transitions. For p < 1 SSL be- 
haves exactly like a regular Neel antiferromagnet, while 
for p > 1 the spin-flop phase is replaced by the umbrella 
phase. The dotted line denotes the special degenerate 
case p = 2. Magnetization curve M(h) are linear until 
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FIG. 5. T = phase diagram of the classical SSL. The solid 
lines correspond to second order phase transitions. The dot- 
ted line mark the special degenerate case J2/J1 = 2. 



saturation, and there is no magnetization plateaus for 
the isotropic Heisenberg SSL. Note that the 1/3 pseudo- 
plateau observed in Ref. 13 was a finite-temperature ef- 
fect. The phase diagram in Fig. 5 has been constructed 
by comparing energies of different spin structures intro- 
duced in this section and selecting the one with the lowest 
energy for a given (p, h). In order to prove that this phase 
diagram is indeed correct, we have to show that there are 
no other spin structures with lower energy. This is done 
in the next section. 



III. EXACT GROUND STATE 

In order to prove that the three spin structures in- 
troduced in the previous section (spin-flop, spin-flip and 
umbrella) are indeed the ground state structures for re- 
spective (p, h) we use the method of decomposing the 
Hamiltonian into overlapping elementary clusters, which 
has been used previously in, e.g., Refs. 1, 13, and 19 



(12) 



A 



For SSL such elementary cluster has the shape of the 
right triangle, highlighted in red in Fig. 1. Each triangle 
includes one spin So at the right angle, two spins (Si 
and S 2 ) at the 45° angles; two Ji bonds, and one J 2 
bond shared by two triangles. The number of triangles 
on the lattice is the same as the number of lattice sites, 
as each triangle includes 3 spins and each spin is a part 
of 3 triangles. The Hamiltonian of a triangle is 



H A = S (S 1 + S 2 ) + |S 1 S 2 -h 



aS + ^^(S 1 + S 2 ) 

where a is an arbitrary real number. It corresponds to 
an arbitrary way in which the term — hSi of the origi- 



nal Hamiltonian can be divided between three different 
triangles which include the site Sj. The triangle Hamilto- 
nian is invariant under three basic symmetry operations: 
simultaneous rotation of all spins around z axis, reflec- 
tion of all spins in xz plane, and interchange Si -O- S 2 , 
where z-axis is parallel to h. These operations generate 
a symmetry group, which also includes such operations 
as reflection in yz plane (or any other plane containing 
z axis) and inversion of the xy components of all spins 
Sf — > -Sf , Sf -> -SV (i = 0, 1, 2). This symmetry leads 
to the degeneracy of most energy levels e(a, So, Si, S 2 ), 
as the only configurations invariant under the symme- 
try group of T-La are collinear ones with So = ±z and 
Si = S 2 = ±z. At the special point p — 2 there is an 
additional symmetry operation So <-» Si (or So S 2 ). 

The triangle Hamiltonian has an a-dependent ground 
state energy eo(a), thus for every possible configuration 
of three unit vectors So, Si, and S 2 the inequality 



6a(S , Si, S 2 , a) > e (a) 



(14) 



holds true. The equality is achieved for a possibly degen- 
erate ground state configuration of the triangle Hamilto- 
nian. There is also an inequality for the lattice energy, 
which holds true for all values of a 



Ne{S,} = J2 6A - J2 e o( a ) = N ^o(a), 



(15) 



A 



A 



or e {Si} > eo(a), where TV is the number of lattice sites, 
and e{Si} is the energy per site. The equality here is 
possible only for a spin structure {Sj} that minimizes 
the energy of each triangle simultaneously. Here and in 
the following we use the word " structure" for spin struc- 
tures {Sj} on the lattice, and the word "configuration" 
for configurations So, Si, and S 2 of the three spins of a 
triangle. According to (15), for every other spin structure 
{S ■} one can write an inequality e {5-} > eo(a) = e {Si}, 
which proves that the structure {Si} is indeed the ground 
state of the lattice Hamiltonian (2), or, in general, one of 
the degenerate ground state structures. In other words, 
in order to find the ground state of the lattice Hamilto- 
nian (2) for a given (p, h), we have to construct a lattice 
spin structure {Si} from the ground state configuration 
(S , Si, S 2 ) of a triangle (let us call it "brick"), or from 
a set of such bricks in case of degeneracy. In case of 
a degenerate ground state, it is important to note that 
any possible lattice structure constructed from bricks is 
a ground state structure, and, vice versa, any possible 
ground state structure can be constructed from bricks 
(which can be seen from the fact that for any ground 
state structure the energy of each triangle is equal to e ). 

The question is whether it is possible to construct a lat- 
tice structure from a given set of bricks, which includes 
all possible realization of the degenerate ground state 
of Ha, spawned by the symmetry group of T-La- The 
problem is not trivial, since each spin is a part of three 
different triangles, so the bricks must match each other 
perfectly. This is obviously not possible for arbitrary a, 
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Energy 



2 + p/2 =F ft 
-2 + p/2 =f ft(l - 2a) 

—p/2 =p fta 
-l/p-p/2-^ 
■ £ + ft fa - ^ ) - ^ 



4p 



P _ 
2 



ft fa 

- a )h 2 



4p 



+ ■ 



Domain 



„>l >h< 2j£±!21 

^<p-l 
l-p< ^ < ! + P 



ft < 



2 
eta' 



TABLE I. Energies e(a) and domains of existence of different 
configurations of a triangle, which are energy extrema of Ha ■ 
We use the definition a' = 1 — a. The umbrella configuration 
is only defined for a = l/(p + 1). 



however, as we are going to see below, such construction 
can indeed be performed for the right choice of a. The 
required values are a = 1/2 for p < 1, and a = l/(p + 1) 
for p > 1, respectively. For these values of a we can 
construct the spin-flop, spin-flip and umbrella structures 
from the ground state configurations of and confirm 
the phase diagram of Fig. 5. Contrary to the case of the 
Ising SSL 19 , the phase diagram of the classical Heisen- 
berg SSL does not have the convexity property, so the 
search for ground state must be performed for arbitrary 
(p, ft), rather than for a finite number of special points. 

For simplicity we assume that ft > 0, so that the field h 
provides a fixed direction z. The rather trivial case ft = 
(solved in Ref. 1 by the same method) has been discussed 
in the previous section. Let us find all possible steady 
states (energy extrema) of the triangle Hamiltonian (13) 
and their energies (listed in Table I). They can be found 
from the system of three vector equations 



_d_ 

OS, 



for the three unit vectors So, Si, S2: 

asT + 2 ~2~ * lSl> 



dn A „ p„ i-\ , „ 

__S +2Si-—h-£ 2 S 2) 



dU 



A 



as 



Si + S2 - ah = 6 So, 



(17) 



(18) 



(19) 



where & are the three real Lagrange multipliers used to 
enforce the conditions S,Sj = 1. They are nonpositive for 
local energy minima, and nonnegative for local maxima. 

First, let us consider collinear (Sj = ±z) solutions of 
Eqs. (17)-(19). Any collinear configuration is a solution. 
The energies of two spin-flip-like configuration (Sj = ±z) 
are tFUpi.i = 2 + p/2 =p ft and do not depend on a. The 
two Neel-like configurations (Si = S2 = —So = ±z) have 



energies ejv ee n,2 = — 2 + p/2 =p ft(l — 2a). Finally, the 
two dimer-like configurations (Si = — S 2 = z, S = ±z) 
have energies CDimeri,2 = —p/2 T ha. Now let us find 
all noncollinear solutions of Eqs. (17) -(19). Subtracting 
first two equations gives 



(20) 



It means that either Si and S 2 are collinear, or that both 
expressions in parentheses are equal to zero. There are 3 
possible cases 

Case 1: 6=6 = -p/2, 



Case 2: Si = -S 2 , 



Case 3: Si = S 2 . 

Below we shall consider all 3 cases in detail. 

Case 1: 6=6 = — p/2- We introduce a new variable 
S = (Si+S 2 )/2, and note that < |S| < 1, but |S | = 1. 
The equations (17)-(19) become 



S +pS = h(l-a)/2, 



2S - ah = 6S - 



Excluding So gives 



S(2 + 6p) = h[a + 6(l-<*)/2] 



(21) 



(22) 



(23) 



There are again two possible cases. 

Case 1.1: 2 + £ p = a + 6(1 - a)/2 = 0. This is only 
possible for a = l/(p + 1), which shows that it is indeed 
the correct value of a for p > 1, since no other value of 
a can produce the umbrella phase. For this value of a 
there is an entire family of the degenerate " umbrella-like" 
solutions of the form 



So+pS 



ph 



2(P+1) 



for ft < h r 



2 (p + l) 2 I p with the energy 



(■Umb 



1 
P 



P 
2 



h 2 p 



4(P + 1) 



2 • 



(24) 



(25) 



For case p = 2, this degeneracy is preserved at the lattice 
level 13 , namely, every configuration with So + Si + S 2 = 
h/3 can form a lattice structure consisting of three types 
of sites (for example, ordered as in Fig. 2, lower right), 
since So, Si, and S 2 can be interchanged freely. In par- 
ticular, the UUD structure can be realized at a single 
point (p, ft) = (2, 3). For p ^ 2 the requirement of bricks 
matching each other partly lifts the degeneracy and leads 
to the umbrella structure, for which Sih — cos is equal 
for all spins. For a = 1/2, p < 1 there is no solutions of 
this kind. 



Case 1.2: S, So and h are collinear (but Si. 2 and h 
are not). This gives up to two Y-likc configurations, with 
S || z and So = ±z. The Y\ configuration with So = z 
has the energy 



EYi (a) 



P 2 



1 - a 



(1 -afh 2 
4^ 



(26) 



and it exists for h < 2(p— 1)/(1— a). The Y2 configuration 
with So = +z has the energy 



1 P u ( 

p 2 V P 



(1 -q) 2 /t 2 
4^ 



(27) 



and it exists for 1 — p < (1 — a)h/2 < 1 + p. 

For a = l/(p + 1), p > 1 the Yi^ solutions are just 
two special cases of the umbrella-like solution introduced 
above. For a — 1/2, p < 1 there is a single Y 2 solution 
with the energy 



1 P 

e Y 2 = r 

p 2 



Hp -i) 

2p 



Ap 



(28) 



Case 2: Si = — S2. The dimer-like collinear solu- 
tions belong to this case, but there are no noncollinear 
solutions. 

Case 3: Si = S2 = S and £1 = £2 = £. The equations 
(17)-(19) take the form 



So 



1 



2S-ah = £oS , 



or, after excluding So, 
£o(l -a) 



a + 



= S 



(29) 



(30) 



(31) 



As before, it can mean either that S is collinear with h 
(which leads to collinear spin- flip and Neel-like solutions), 
or that both expressions in the square brackets are equal 
to zero, which leads to 



£0 = -: 



2a 



^ 2 a 



and 



S + Sr 



a 



I- a 2 

which is the spin-flop-like solution with the energy 



eFiop(a) = - 



a(l 



1 



a 



(32) 



(33) 



(34) 



It exists for h < 2/(a(l - a)). For a = l/(p + l),p > 1 
this is again a special case of the umbrella-like solution. 
For a = 1/2, p < 1 this is the spin- flop solution (S + S = 
h/4) with the energy e F lop = -2+p/2 — h 2 /l6. The spin- 
flop-like bricks only match if Sh = Soh, i.e. that all spins 




FIG. 6. (Color online) Magnetization curves M(h) for the 
classical SSL. p = 0.5, no anisotropy: exact result (solid black 
line), numerical data (solid circles), p — 2.5, no anisotropy: 
exact result (dotted black line), numerical data (empty cir- 
cles), p = 2.5, with anisotropy: f3 = 0.1 (blue dashed curve), 
P — 0.3 (red dash-dot curve). All numerical calculations have 
been performed on the 12 x 12 lattice with periodic boundary 
conditions. 



are tilted by the same angle 8' relative to the xy-plane. 
This condition gives a = 1/2 as the only value of a for 
which the construction of the spin-flop lattice structure 
is possible. 

Some of the configurations corresponding to the energy 
extrema of Ha do not match and thus cannot form a spin 
structure on the lattice. The ones that do are: spin-flop 
and Neel (for a = 1/2), umbrella- like (for a — l/(p+l)), 
dimer (for a = 0), and spin-flip (for any a). The direct 
comparison of their energies (listed in Table I) for respec- 
tive values of a shows that spin-flop, umbrella-like and 
spin-flip solutions indeed minimize the triangle Hamil- 
tonian (13) in the three respective regions of Fig. 5. 
The Neel structure is only realized for h = 0, p < 1. 
The dimer structure does not exist for isotropic classical 
Heisenberg SSL. This confirms the phase diagram shown 
in Fig. 5. 



IV. NUMERICAL SIMULATIONS AND THE 
EFFECT OF ANISOTROPY 



In order to give an independent check of our exact 
results we performed a series of numerical simulations, 
calculating the ground state of lattice Hamiltonian H for 
different values of p, h. We have also examined the effect 
of uniaxial anisotropy by adding the term 



1 - (si 



(35) 



7 



to the Hamiltonian %, or, equivalently, adding the term 

;i - {sif) (36) 



to where (3 = B/J\ is the anisotropy constant. 

The minimum of % was found by a discrete micro- 
magnetic simulation with only the Gilbert damping term 
included. We considered the system of equations 



dt 



= -A 



&H_ 



(37) 



where A > is the damping parameter, for each spin 
of the lattice; and solved it using the first-order Runge- 
Kutta method in spherical coordinates with our own 
computer code. This method decreases the energy of the 
system on each step (provided that A is small enough) 
eventually finding a (local) minimum. It can be viewed 
as a variation of the steepest descend method. Unlike 
previously used Monte Carlo methods, our method looks 
for the ground state of the system avoiding any finite- 
temperature effects. We have used cells of different size 
(6 x 6, 12 x 12, and 24 x 24) with periodic boundary 
condition. In order to minimize the probability of find- 
ing a local energy minimum instead of the ground state, 
each simulations was performed 50 times with different 
random initial conditions, and the result with the lowest 
energy was chosen. 

A number of simulations for different values of p, h 
have been performed. Our results fully confirm the phase 
diagram of Fig. 5, in particular, the second order phase 
transitions and the additional degeneracy for p — 2 are 
clearly seen in our calculations. We have also applied 
the same method to the triangle Hamiltonian and 
checked numerically the validity of the results of the pre- 
vious section. Two typical magnetization curves M(h) 
are presented in Fig. 6 for p = 0.5 and p = 2.5. The nu- 
merical results are plotted as circles in this figure, while 
the analytical formulae (4)-(10) are presented as lines. 
For p = 0.5 the system goes through the sequence of the 
Neel-spin-flop-spin-flip structures, all of them being pe- 
riodic with 4 atoms per unit cell, exactly as predicted by 
our analytical treatment above. The M(h) curves in fact 
do not depend on p at all as long as p < 1. For p > 1 
the ground state is the umbrella structure, which is in 
general incommensurate with the lattice. The numerical 
calculations with a finite lattice size cannot reproduce 
this phase exactly, of course, but as long as the lattice 
size is 6 x 6 or larger, the resulting umbrella-like phase 
has magnetization and energy very close to the values 
given by Eqs. (10), (11). Even for the 6x6 lattice the 
relative difference between theoretical and calculated M 
is only about 9 x 10~ 4 , a difference too small to be seen in 
the scale of Fig. 6. Examples of such umbrella structures 
with periodicity forced by the lattice size are presented 
in Fig. 3 and Fig. 4 for (p,h) = (1.5,3) and (1.2,3) 
respectively. 



If the magnetic anisotropy /3 is switched on, the single 
UUD point (p, h) = (2, 3) expands into a finite region of 
the UUD phase. The M(h) curves for p ^ 2 first get a 
noticeable kink at around M = 1/3 (Fig. 6, dashed blue 
curve), which eventually turns into the M = 1/3 UUD 
plateau when the anisotropy constant is increased (Fig. 
6, dashed-dot red curve). For p = 2 the M = 1/3 step 
appears for any finite value of f3. Our M(h) curves are 
very similar to the M(h) curves of Ref. 14, which were 
obtained using Monte Carlo method for p = 2. However, 
we observed the onset of the 1/3 plateau at a finite value 
of f3, which was not seen in Ref. 14 due to the choice of 
a non-arbitrary point p = 2. No plateaus other than 1/3 
appear in our calculations. 

It is interesting to compare the different physics of Ising 
and classical Heisenbcrg SSL. The Ising SSL allows for 
collinear structures only, and its zero-temperature M(h) 
curves consist of vertical and horizontal segments only 
(plateaus and first order phase transition). In contrast, 
magnetization curves for isotropic Heisenberg SSL are 
linear until saturation thanks to noncollinear spin-flop 
and umbrella structures. The anisotropic Heisenberg SSL 
combines features of both Ising and Heisenberg models. 
Its magnetization curves M{h) can include both hori- 
zontal steps (corresponding to collinear spin structures) 
and inclined regions (corresponding to noncollinear struc- 
tures). While the Ising approach might be sufficient for 
T111B4, both types of regions are clearly seen experimen- 
tally in H0B4. In fact, the experimental M(h) curve for 
H0B4 (Ref. 12) is remarkably similar to our M(h) curve 
for p = 2.5, (3 = 0.3 (dashed-dot red curve in Fig. 6). 
The next logical step in the study of classical Heisenberg 
SSL would be determining of the phase diagram of the 
classical Heisenberg SSL with easy-axis anisotropy. This 
problem is, however, beyond the scope of the present pa- 
per and will be addressed in our future research. 



V. CONCLUSION 

We have determined analytically the exact ground 
state of the classical Heisenberg SSL in the external field 
h and presented the phase diagram of this model in the 
(p, h) plane. The phase diagram includes the spin-flop 
phase for p < 1, the umbrella phase for p > 1, and the 
spin-flip phase for sufficiently large h. The phase tran- 
sitions between these three phases are of the second or- 
der. The zero-temperature magnetization curves M(h) 
are linear until saturation with no features. For p = 2 
there is an additional degeneracy and an infinite number 
of spin structures which share the energy and magnetiza- 
tion with the umbrella one. In particular, UUD structure 
can be realized at a single point (p, h) = (2,3). 

The effect of the easy-axis uniaxial anisotropy on M(h) 
curves has been examined by numerical micromagnetic 
simulation. The anisotropy leads to the onset of the 
M = 1/3 UUD plateau at a certain finite value of the 
anisotropy constant /3 for p ^ 2 (and at j3 = for p = 2). 
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Our results demonstrate the existence of both collinear 
(steps of the M(h) curve) and noncollinear (inclined parts 
of the M{h) curve) spin structures for anisotropic Heisen- 
berg SSL. The results explain the magnetization curve of 
H0B4, while the M(h) curves of TmB4 are more Ising-likc 
in nature, and cannot be explained without introducing 



additional interactions into the model. 
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